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A systematic pointing error model is used to caiibrate antennas in the Deep 
Space Network. This article describes and analyzes the least-squares problem and 
the solution methods used to determine the modeTs parameters. Specifically studied 
are the rank-degeneracy problems resulting from beam-pointing error measurement 
sets that incorporate inadequate sky coverage. A least-squares parameter subset 
selection method is described and its applicability to the systematic error modeling 
process is demonstrated on a Voyager 2 measurement distribution. 


I. Introduction 

A pointing error model is used in the Deep Space Net- 
work’s (DSN’s) antenna-calibration process. With the ex- 
ception of environmental effects, the major sources of er- 
rors in an antenna-pointing system are systematic and 
repetitive and therefore can be closely modeled. Exam- 
ples of parameters in the model are residual errors in the 
geometric alignment of the mount axes and fixed-angle en- 
coder offsets. Data collected from spacecraft and radio star 
observations are used to determine the parameters in the 
model and are then entered into the pointing system to 
accurately point the antenna. The origins of the pointing 
error modeling approach for radio- frequency (RF) anten- 
nas can be found in [1,2] while its development within the 
DSN is discussed in [3]. 1 


1 R. L. Riggs, “Antenna Pointing Angie Corrections,” DSN Antenna 
Seminar, Videotapes 49-54, Jet Propulsion Laboratory, Pasadena, 
California, May 1986. 


The complete pointing error model is the sum of its 
separate error components. Table 1 shows individual error 
sources and the elevation/cross-elevation (or declination/ 
cross-declination, depending on antenna mount) regressor 
variables used to estimate parameters. See [1,2,4] for a 
more thorough discussion of these parameters. Currently, 
this entire model is set in motion in the antenna-pointing 
system by entering parameter values manually. The DSN 
70- m antennas track targets in both the computer com- 
mand and precision modes of operation, each defined by a 
set of axis position transducers. (See [5] for a discussion 
of the axis servos and controllers.) The 34-m antennas 
employ only the computer command mode. In [3,6] rec- 
ommended model parameter sets are given that apply to 
each tracking mode of these antennas; they are also re- 
peated in Table 2. As can be seen, nine error parameters 
are used to estimate in the precision mode and eight in 
the computer command mode. In practice, the model pa- 
rameters are determined by performing a least-squares fit 
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on the pointing offset data collected from the spacecraft 
and/or from radio star observations. In this article, refer- 
ence will only be made to the particular combinations of 
parameters in Table 2. 

This article explores the numerical properties of the sys- 
tematic error modeling process. Specifically, the analysis 
focuses on the numerical properties of the matrix formed 
by the pointing model regressor variables evaluated over 
the beam-pointing error data sets. These measurement 
sets may not cover enough points in the sky to accurately 
estimate all of the parameters. This is due to the finite 
number of targets and to other practical operational con- 
straints, such as lack of antenna time. On the other hand, 
the objective of particular calibrations may be to opti- 
mize pointing in a particular region of the sky, such as 
along a constant declination. In practice, however, the lim- 
ited measurement sets lead to rank deficiency in the least- 
squares measurement distribution matrix. This study of 
the problem will lead to a more objective approach to pa- 
rameter selection and parameter estimate interpretation. 
In addition, the analytical techniques provided here may 
be used to predict which directions in the sky will yield 
optimal estimation. 

The remainder of this article will formulate the system- 
atic error parameter estimation problem and then estab- 
lish a hypothetical performance index for matrix condi- 
tioning. In addition, the numerical tools presented will be 
used to analyze practical sky distributions in the context 
of the least-squares approximation and the current solu- 
tion method will be reviewed. The article concludes with 
a proposed algorithm for parameter selection. 


II. Model Generation 

In order to accurately point the antenna, pointing er- 
ror correction models must be generated from radio star 
or spacecraft pointing offset data. This section deals with 
the model-fitting process, which uses the least-squares al- 
gorithm and assumes that the measurement data sets are 
accurate. At this time, the estimation process does not 
deal with uncertainties in the conical scan pointing offsets 
and radio star boresights except in human filtering of very 
large nonrepeatable and unexplainable offsets. 

A. Least-Squares Problem Formulation 

The parameter vector p of the systematic pointing error 
model is determined by performing a linear least-squares 
fit on the offset data. The estimation problem is formu- 
lated from m observations as 


" Sxel i ' 


Axeli(eli, azi) 

Sxelrn 


Axel m (^el m , fl£ m ) 

del 1 


Aeli(el u azi) 

- - 


. Ae/ m (e/ m , nz m ) 


or 

y = Ap (2) 

where the offset vector y is 2m x 1, the measurement dis- 
tribution matrix A is 2m x n, arid the parameter vector p 
is n x 1. As can be seen, equations representing both the 
cross-elevation and elevation error functions are obtained 
for each single observation point in the sky. Let the least- 
squares estimator be p and satisfy the following matrix 
equation 

y = Ap (3) 

where the vector y contains the estimated (or fitted) values 
to the cross-elevation and elevation offsets of Eq. (1). The 
difference between the individual elements, or residuals, is 
defined as 

r i = Vi ~ Vi (4) 

The method of least squares chooses the parameter esti- 
mate p, such that the following quantity is minimized 

2m 

X>? (5) 

i—l 

The estimate satisfies the following matrix equation [7,9] 

p = (A t A)~ 1 A t y (6) 

where A* is the transpose of A. Caution must be given 
to least-squares problems in which the regressor variables, 
or basis terms of A, are not truly independent. In such 
cases the measurement distribution matrix A may be close 
to, or is, rank deficient. If A is rank deficient, then there 
are an infinite number of solutions to the least-squares 
problem and no conclusion can be drawn as to the role of 
the individual regressor variables [9]. 
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During the systematic pointing error estimation pro- 
cess, limited data sets and inherent correlations in the 
pointing error model have led to rank deficiency and its 
associated problems. This situation was discussed in [8] 
where the condition of empirically correlated regressor 
variables was termed “multicolinearity.” It was pointed out 
in [8] that regressor variables of the model are not truly 
independent. However, this is not accurate when a proper 
combination of parameters is selected, as recommended in 
Table 2. The degree of linear independence in the columns 
of matrix A for various antenna configurations is strongly 
dependent on the distribution of the observation points 
over the sky. This situation and its effect on the pointing 
error estimation is discussed below. 


and 

<7i > Oo > ■ • * > &n > o (9) 

The theorem is taken from [10] and the more gen- 
eral SVD is proven in [9]. The numbers oq , (To, . . . , <r n , 
which are unique, are called the singular values of A. The 
columns [u 1 ,u 2 , . . . ,u m ] of U are called the left singular- 
vectors of A, and the columns [vi,v 2 , v u ] of V are 
called the right singular vectors of A. SVD is extremely 
useful in analyzing numerical rank deficiency because the 
singular values indicate how near A is to a matrix of lower 
rank. The matrix A has rank r if and only if 


B. All-Sky Model Analysis 

An analytical approach was taken to obtain a perfor- 
mance index for the numerical conditioning of the sys- 
tematic error least-squares problem and to compare it with 
results from practical measurement sets. One such perfor- 
mance index can be determined by examining a hypothet- 
ical all-sky uniform distribution of pointing offset data. 
These measurement points are used to generate measure- 
ment distribution matrices for different combinations of 
parameters. Intuitively, it would make sense to obtain 
pointing error offsets uniformly throughout the field of 
view of the antenna and conclude that this is the opti- 
mal distribution for input into the parameter estimation 
problem. However, observing the basis terms of the point- 
ing model given in Table 2, it can be seen that not all 
terms are simultaneously functions of both azimuth and 
elevation. This condition will tend to result in redundant 
column elements of A; thus, optimal matrix conditioning 
will most likely not be obtained with the all-sky distribu- 
tion. However, as will be shown, all-sky matrices do have 
acceptable conditioning and can be used for a suitable per- 
formance index. Singular value decomposition (SVD) was 
used to analyze the linear independence of the columns of 
A and is defined in the following theorem. 


<?r > 0 = cr r + 1 (10) 

Mathematically speaking, the smallest singular value 
of A is the 2-norm distance of A to the set of all rank- 
deficient matrices [9]. The ratio of the largest to small- 
est singular value is termed the condition number of A. 
This number quantifies the sensitivity of the least-squares 
solution p of Eq. (2). Large condition numbers indicate 
that relatively small changes in A or the offset vector y 
can induce large changes in the computed least-squares 
solution p. This is undesirable since parameter vector 
estimates computed from such ill-conditioned measure- 
ment distribution matrices can lead to erroneous pointing 
offset corrections that will be applied in future antenna 
tracks. The pointing model corresponding to the DSS 14 
70- m antenna (i.e., latitude = 35.426) was used to gener- 
ate the full A matrix of relevant error parameters. The 
uniform distribution consists of 9-deg increments in eleva- 
tion and 20-deg increments in azimuth. This full measure- 
ment set is illustrated in Fig. 1. The singular values of 
the A matrices corresponding to the precision mode and 
computer command mode of operation are presented in 
Table 3. 


Theorem 1. Let A be a real m x n matrix with m > n. 
Then there is an orthogonal matrix U = [ui , . . . , u m ] of 
order m and an orthogonal matrix V = [vj , . . . , v n ] of or- 
der n such that 


U l AV 


£ 

0 


where 


£ = diag(<ri,...,<r n ) 


(7) 


( 8 ) 


As can be seen, the numerical conditioning for the least- 
squares problem resulting from this hypothetical all-sky 
distribution is well behaved. Both the precision mode and 
computer command mode parameter sets yield A matrices 
with reasonably nonzero singular values and small condi- 
tion numbers, implying full-column rank. Another quan- 
tity commonly considered in least-squares analysis is the 
correlation matrix derived from (A 1 A ) -1 of Eq. (6), which 
is numerically shown in Table 4. The matrix (A* A ) -1 is 
an estimate of the covariance matrix for the solution vector 
of the least-squares problem. Values near one in the cor- 
relation matrix indicate high pairwise correlation between 
the estimated parameters. 
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It is evident that such high correlation is implied be- 
tween the first three parameters (az collimation, fixed az 
encoder offset, and az/el axis skew) of the computer com- 
mand mode set. Evaluation of the basis terms correspond- 
ing to these parameters in the limited elevation range of 0 
to 90 deg results in the pairwise correlation and cannot be 
avoided regardless of the azimuth distribution. The impli- 
cation of inherent correlation to parameter estimate sta- 
bility was investigated through Monte Carlo simulations. 
The empirical estimation covariance matrix was computed 
and found to be in very close agreement with the theoret- 
ical covariance matrix computed from (A 1 A)" 1 , thus il- 
lustrating that stable parameter estimates will result from 
an all-sky distribution. It was noted in the simulations 
that individual estimates of the first three elements al- 
ways varied in the same direction of magnitude, but that 
differences never exceeded the bounds predicted in the the- 
oretical standard-deviation vector given by 

n 


where the v\ and <r,- are defined in Theorem 1. The above 
equation is obtained by solving for A in Eq. (7), substi- 
tuting it into (A t A)~ 1 , and then taking the square root 
of the diagonal of the resulting matrix. 

The numerical conditioning of the least-squares estima- 
tion of antenna precision and computer mode systematic 
error parameter sets was evaluated above for an all-sky 
distribution. The resulting measurement distribution ma- 
trix for each mode of operation was found to have full 
rank, thus ensuring unique least-squares solutions for p. 
Also, the large values in the correlation matrix were not 
seen to degrade the stability of repeated parameter esti- 
mates. The linear dependence of the parameters implied 
by the correlation matrix is due to their mathematical defi- 
nitions and selecting them simultaneously will not degrade 
the estimate of the measurement vector y. Such rich off- 
set distributions can never be obtained in practice, thus 
it is inevitable that poorer matrix conditioning will lead 
to least-squares estimates of poorer quality. As shown by 
this analysis, the singular values and condition number of 
the distribution matrix A are key parameters in evaluating 
ill-conditioned least-squares problems. 

C. Reduced and Sparse Data Sets 

Current practices dictate that systematic error models 
be generated from antenna-pointing error-correction data 
taken from as much of the sky as possible or from an area 
defined by one or two declination angles. The first is used 


to generate all-sky pointing models, while the second com- 
putes model parameters applicable only in limited direc- 
tions of the sky. Both situations typically diverge from the 
hypothetical all-sky example since the basis terms of the 
pointing model are evaluated in fewer, and perhaps more 
redundant, directions. Their effect on the least-squares 
estimation process will be illustrated with examples. 

Figure 2 shows the sky trajectory for the Voyager 2 
spacecraft. Conical-scan offset data collected at a decli- 
nation of —22.5 deg clearly represent only a small por- 
tion of the total sky measurement space. Tables 5 and 
6 show the singular values, condition numbers, and theo- 
retical standard deviations in millidegrees (mdeg) for the 
least-squares estimate using the A matrices generated for 
precision and computer command operation. As implied 
in the tables, matrix condition deteriorates in both param- 
eter sets because of reduced measurement space. The the- 
oretical standard deviations of the all-sky parameter sets 
are shown in mdeg in Table 7. Comparison with those 
of Table 6 illustrate the degradation of the least-squares 
parameter estimation. In [8], least-squares parameter fits 
were done on Voyager 1 conical-scan data obtained from 
the DSS 14 64-m antenna. The results were parameters 
that were too large in magnitude to be realistic or practi- 
cal and that were unstable on a day-to-day basis. It has 
been shown through SVD analysis that such ill condition- 
ing of the systematic error least-squares problem can, in 
general, be inferred a priori for any constant declination 
measurement set. 

Figure 3 shows sky distribution for a radio source bore- 
sight offset file taken at the DSS 13 26-m antenna. The dis- 
tribution is typical of data gathered during planetary radio 
astronomy experiments — here, for four radio sources. The 
pointing model regressor values were once again evaluated 
at the source coordinates and results of the SVD analysis 
are shown in Tables 8 and 9. Condition numbers for pre- 
cision and computer command mode parameter sets are 
comparably small in magnitude to those from an all-sky 
distribution. The smallest singular values are also reason- 
ably nonzero. Only minimal estimate degradation is pre- 
dicted by the increase in theoretical standard deviations. 
Furthermore, the magnitude of this uncertainty is still rea- 
sonably small in the context of parameter estimates, which 
are usually in the tens of mdeg. This example illustrates 
that rank deficient measurement distribution matrices can 
be avoided by using recommended parameter sets and by 
evaluating the regressor variables with adequate sky dis- 
tribution of pointing offsets. Concluding that such a mea- 
surement set is adequate for the least-squares model fitting 
is essentially putting emphasis on the norm of the result- 
ing solution vector p instead of minimizing the norm of the 
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residual vector r of Eq. (4). This approach appears to be 
the most logical given that the measurement uncertainties 
are not modeled. It has also been shown that, for radio 
source pointing calibrations, this matrix condition analy- 
sis can be done during pretrack activities, thus influencing 
the scheduling of calibrators. 

It must be stressed that these results hold only for the 
sets of parameters recommended in Table 2. Different com- 
binations of 21 error coefficients in the current pointing 
model will yield different, and in some cases disastrous, 
numerical properties of the matrices involved in the com- 
putation of the least-squares solution. 


III. Solution Methods 

A. Parameter Selection 

The two goals of the modeling process are to quan- 
tify contributors to the antenna’s systematic pointing error 
so that pointing can be corrected and so that knowledge 
of the antenna’s mechanical and structural characteristics 
can be acquired. To achieve both of these objectives si- 
multaneously, identical parameter vectors must be chosen 
for estimation on a consistent basis. These parameters for 
the 70-m and 34-m antennas in the applicable mode of op- 
eration have been given in Table 2. Subsets of these vec- 
tors should be chosen either when parameter values are 
physically known a priori or when they are consistently 
estimated with small magnitudes. In practice, however, 
the goal of correct pointing can be achieved without ac- 
curate knowledge of actual antenna error characteristics. 
Optimization may be based on any random set of parame- 
ters that minimizes the sum of the squares of the residuals 
given in Eq. (5) without regard for physical interpretation. 

Regardless of the estimation philosophy practiced, 
problems always arise when building models for partic- 
ular regions of the sky — for example, along a line or band 
of constant declination for one or more sources — or for 
a particular spacecraft. The rank deficiency that plagues 
least-squares problems in these cases generates uncertainty 
in parameter selection and interpretation. However, such 
models for locally optimized pointing are needed for criti- 
cal spacecraft and holography tracks and for those tracks 
of single sources known as strong, reliable antenna calibra- 
tors. The current least-squares solution method described 
next uses the SVD to accommodate ill-conditioned mea- 
surement distribution matrices. 

B. Singular Value Decomposition 

The SVD subroutines in the systematic error modeling 
software that were used to solve the least-squares problem 


were taken from [11]. A key feature of the SVD method 
is its ability to handle rank deficiency. Ill-conditioned A 
matrices result in the rank (A) = r being less than the 
parameter dimension n. This results in a rank-deficient 
least- squares problem that has an infinite number of so- 
lutions, for if the vector p is a minimizer and the vector 
z € null( A), then p + z is also a minimizer. The SVD 
method is useful in such situations since it is a reveal- 
ing and complete orthogonal decomposition. The routines 
from [11] basically implement the following theorem taken 
from [9], given here without proof. 

Theorem 2. Suppose U t AV = X is the SVD of 
A € ® mxn with r = rank(A). If U = [uj , . . . , u m ] and 
V — [vi, . . . , v n ] are column partitionings and y £ 3? 111 
then 

PLS =£^ V| ( 12 ) 
<=i ai 

minimizes || Ap — y 1 1 2 and has the smallest 2-norm of all 
minimizers. Moreover 


II Ap LS - y H2 = ( u ‘y ) 2 ( 13 ) 

t = r + l 


Note that if r < ?i, this corresponds to simply adding a 
zero multiple to the solution vector pls rather than adding 
random large-valued multiples produced by the near-zero 
singular values. This may reduce uncertainty in the esti- 
mated coefficients, as in Eq. (11), but increases the resid- 
ual norm, as in the increased summation index of Eq. (13). 
This point was touched on earlier. In practice, one must 
still come up with a numerical estimate r of r. The sys- 
tematic error modeling software estimates the numerical 
rank r of A as 


' (Tr > 6 > <Ty + 1 > • • ■ <r„ ( 14) 


where the tolerance 8 is chosen to be < 7 2 7 scaled by a 
machine-precision dependent factor. The selection of r 
results in 


Pr = 


1 = 1 



(15) 


as an approximation to pls If oy >> 8, then p;. is a 
very close approximation to the true minimizer pls since 
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A can be unambiguously regarded as a matrix with rank 
r [9]. When [<?i , . . . , <7„] do not clearly split into small 
and large values, rank determination may be somewhat 
arbitrary. 

C. A New Algorithm for Parameter Selection 

As has been shown, the SVD solution currently alle- 
viates the rank-deficiency problems associated with lim- 
ited pointing ofTset distributions. This means filtering out 
small singular values of A and replacing them with the 
matrix A* defined as 


r 

Ar = 53<7«UiV* (16) 

1=1 

where r is the numerically determined estimate of the rank 
of A. As discussed in [9] such a cutofT makes sense when 
the measurement distribution matrix is derived from noisy 
data. However, in this case, A is being evaluated using 
accurate ephemeris from observed targets, as in Eq. (1). 
In other applications, rank deficiency is an indication of 
redundancy among factors that comprise the model. As 
has been shown in previous sections, redundancy among 
systematic error regressor variables occurs in the estima- 
tion process only when dealing with limited and reduced 
pointing measurement sets. In these cases, the system- 
atic error predictor A^p* used in subsequent tracking will 
involve all n redundant factors that may have been cho- 
sen as a result of random parameter selection. Although 
such solutions may correct future pointing, parameter esti- 
mates can obscure physical interpretation of true antenna 
mechanical characteristics. In such instances, it is argued 
that the least-squares solution vector should contain at 
most r nonzero systematic error parameters, which in turn 
dictate which columns of A will be used in approximating 
the observation vector y. The problem of choosing the ap- 
propriate columns of the measurement distribution matrix 
is termed subset selection. The SVD-based subset selec- 
tion procedure that has been chosen for this least-squares 
application is summarized below: 

(1) Compute the SVD A = UEV* and use it to deter- 
mine a rank estimate r. 

(2) Calculate a permutation matrix P such that the 
columns of the matrix Bj € 3? mxr in AP = [Bj B 2 ] 
are “sufficiently independent.” 


criteria instead of those used in Eq. (14). The new criteria 
are based on the matrix condition number and the mag- 
nitude of the theoretical standard-deviation vector given 
by Eq. (11). Given r, the first r columns of permuta- 
tion matrix P give the column indices of A for use in the 
least-squares estimation. These are equivalent to the pa- 
rameters from which the model is selected. A thorough 
discussion of the various approaches to this problem can 
be found in [9,10]. Below is a summary of the algorithm to 
compute P that was chosen and implemented in the sys- 
tematic error modeling software. It is based on both the 
SVD and on QR factorization with the column-pivoting 
algorithm. For A £ dl mxn , QR factorization with column 
pivoting from [9] produces AP = QR where 

(17) 


where r is the rank (A), Q is orthogonal, Rjj is upper tri- 
angular and nonsingular, and P is a permutation matrix. 
This factorization implies that the first r columns of Q 
form an orthonormal basis for range (A). It is the desired 
result since the measurement vector y in the least-squares 
problem may be approximated by the first r columns of 
the matrix AP, which is just Bi of the second step above. 
This is equivalent to choosing the first r parameters of Pp 
for estimation, which is equal to the vector z in step three 
above. As in a previous section where the case for cutting 
of! singular values in the SVD method w*as presented, re- 
ducing the order of the parameter solution vector will also 
increase the residual norm. 

Unfortunately, QR factorization with column pivoting 
alone is not a totally robust method for computing the 
permutation matrix P [9]. The preferred algorithm im- 
plemented in the software that uses both SVD and QR 
factorization is presented in the revised steps below: 

(1) Compute the SVD A = USV 1 and use it to deter- 
mine a rank estimate r. Save the matrix V. 

(2) Apply the QR factorization with column pivoting 
to the subset of V 4 : Q 4 V(:, 1 : r)*P = [R n i£ 12 ] 
and set AP = [Bj B 2 ] with B 1 £ 3P nXr and B 2 £ 

Sj^mx(n-r) 

(3) Determine z £ which minimizes || Biz — y || 2 . 


R: 


Rn 

0 


Rl2 

0 

n — r 


(3) Predict y with the vector Ap dU b where p s „b = [z O] 4 
and z £ 5f r minimizes || B i z — y 1 1 2 - 

Using systematic error modeling, the rank determina- 
tion in the first step can be chosen with more heuristic 


The main contribution of this algorithm is facilitating 
parameter selection for reduced and constant declination- 
pointing measurement sets. In the latter case, its appli- 
cation will ensure consistent parameter selection for par- 
ticular radio sources and spacecraft tracks. This subset 
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selection procedure essentially eliminates parameters that 
the algorithm has deemed unobservable in the given mea- 
surement distribution. The next step is to decide how 
to deal with these excluded parameters. One approach 
is to simply ignore them and proceed as usual with the 
least-squares estimation with the reduced vector z, as de- 
termined above. 

A different approach, when possible, is to use physically 
known or accurate a priori estimates for unobservable pa- 
rameters and subtract their contributions from the mea- 
surement vector y before estimating z. Such an option 
is available in the current software. The resultant solu- 
tion vector should be more consistent with all-sky models. 
Finally, it should be noted that the tools presented here 
can be used for the opposite effect (e.g., predicting matrix 
condition and rank or for least-squares estimate accuracy) 
when the measurement vector y is augmented with point- 
ing offsets taken in new directions. 

To illustrate, this algorithm is applied to the A matrix 
which resulted from the Voyager 2 trajectory, as shown 
in Fig. 2. Referring to Table 5, one can base the rank 
determination of the A matrix on the smallest of the sin- 
gular values. For example, choosing 0.1 as a singular value 
cutoff results in precision and computer command mode 
parameter selections and matrix conditions that are sum- 
marized in Tables 10 and 11. Eliminating parameters 11 
and 21 from the precision mode and 1 from the computer 
command mode results in reduced matrix condition and 
smaller estimation standard deviations for some elements 
of the solution vector. In practice, the actual systematic 
error estimated values are generally less than 100 mdeg. 
Thus, estimation accuracy for some of the remaining pa- 
rameters in Table 10 will be a certain percentage of the 
estimated values. 

Depending on the antenna’s frequency band, this may 
or may not meet the pointing requirements. (A detailed 
description of errors will not be covered here.) Estima- 
tion errors will always be larger in practice because of 
uncertainties in the measurement vector y, so one may 
decide to increase the singular value cutoff and apply the 
subset selection algorithm. Using cutoffs 1.0 for the preci- 
sion mode and 0.2 for the computer command mode yields 
the results summarized in Tables 12 and 13. To achieve 
accuracy comparable to the hypothetical all-sky models, 
the parameters to be excluded are 1,7,11, and 21 from 
the precision-mode set and 1 and 7 from the computer 
command-mode set. It is advised that whenever the fixed 
angular encoder error parameters (for example, 7 and 21) 
are excluded in the subset selection procedure, their values 
should be determined directly from the pointing offset data 


and contributions to y should be removed before making 
an estimation. 

IV. Summary 

This article lias described and analyzed the least- 
squares problem encountered in the DSN systematic point- 
ing error modeling process. Specifically investigated is 
the relationship between rank degeneracy of measurement 
distribution matrices and limited-sky distributions of the 
pointing error offsets. Using a hypothetical all-sky per- 
formance index and an SVD analysis, it is shown that an 
acceptable matrix condition of the least-squares problem 
can be obtained by evaluating the pointing model regressor 
variables with adequate sky distributions of the pointing 
measurements. In addition to matrix condition, the the- 
oretical standard deviations of the least-squares estimate 
are used to evaluate accuracy. It is shown through an ex- 
ample that redundancy among the systematic error model 
regressor variables occurs when dealing with limited and 
sparse data sets. In practice, rank-degenerate mat rices are 
encountered when building models for particular regions of 
the sky, such as along a band of constant declination. 

The key feature of the analysis presented is its predic- 
tive capability. Matrix condition and least-squares esti- 
mate accuracy based on measurement distribution may be 
predicted before actual pointing calibration activities com- 
mence. The current least-squares solution method based 
on singular- value decomposition is also presented. This 
method can handle ill-conditioned measurement distribu- 
tion matrices encountered in the model-building process. 
For limited measurement sets, it was argued that it may 
be preferred to estimate only observable parameters. Sys- 
tematically eliminating redundant parameters will facili- 
tate the parameter selection process and make it consis- 
tent. A recommended subset selection algorithm based 
on singular- value decomposition and QR factorization is 
illustrated with a Voyager 2 measurement set. 

V. Future Work 

This article has presented an analytical approach using 
mathematical tools to answer fundamental numerical ques- 
tions arising from the systematic error-modeling process. 
Such general but consistent procedures are needed in the 
modeling process because of the many antenna-specific me- 
chanical and other practical considerations encompassed 
by the problem. Once past this juncture, one may begin 
to address t he deficiencies and look for possible refinements 
in the estimation process. The most obvious is that of re- 
cursive estimation. Methods must be devised to handle 
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delta sets spanning many weeks or years and incorporat- 
ing many a priori models and model uncertainties into the 
estimation algorithms. If and when uncertainties in point- 
ing measurements can be accurately modeled, including 
those from natural or manmade sources as well as from 
antenna-system imperfections, then the algorithm should 
also be modified to allow for weighted observations. 


All these enhancements must be worked into the 
existing modeling software. This package should also en- 
able the functional form of the model to change relatively 
often. This will allow for the addition of newly discovered 
error terms and for enhancements to accommodate new 
antenna architectures such as the DSS 13 beam-waveguide 
antenna. 


Acknowledgments 


The author would like to thank R. Riggs for instruction about the DSN sys- 
tematic pointing error modeling process. The modeling software developed and 
maintained by Riggs for the Goldstone tracking station was used in this study. 
Chris Sanelli is acknowledged for implementing and testing the parameter selection 
algorithm. 


References 


[1] P. Stumpff, “Astronomical Pointing Theory for Radio Telescopes,” Klein- 
Heibacher Benchte, vol. 15, Fornmoldeteshnischon Zentralamt, Darmstadt, West 
Germany, pp. 431-437, 1972. 

[2] M. L. Meeks, H. A. Ball, and A. B. Hull, “The Pointing Calibration of the 
Haystack Antenna,” IEEE Transactions on Antennas and Propagation , 
vol. AP-16, no. 6, pp. 746-751, November 1968. 

[3] R. Stevens, R. L. Riggs, and B. Wood, “Pointing Calibration of the MKIVA DSN 
Antennas for Voyager 2 Uranus Encounter Operations Support,” TDA Progress 
Report 42-87, vol. July-September 1986, Jet Propulsion Laboratory, Pasadena, 
California, pp. 206-239, November 15, 1986. 

[4] C. N. Guiar, F. L. Lansing, and R. Riggs, “Antenna Pointing Systematic Error 
Model Derivations,” TDA Progress Report 42-88, vol. October-December 1986, 
Jet Propulsion Laboratory, Pasadena, California, pp. 36-46, February 15, 1987. 

[5] R. E. Hill, “A Modern Control Theory Based Algorithm for Control of the 
NASA/JPL 70-Meter Antenna Axis Servos,” TDA Progress Report. 42-91 , 
vol. July— September 1987, Jet Propulsion Laboratory, Pasadena, California, 
pp. 285-294, November 15, 1987. 

[6] C. N. Guiar, R. L. Riggs, R. Stevens, and M. Wert, “DSS 14 Antenna Cal- 
ibrations for GSSR/VLA Saturn Radar Experiments,” TDA Progress Report 
42-93 , vol. January-March 1988, Jet Propulsion Laboratory, Pasadena, Califor- 
nia, pp. 309-337, May 15, 1988. 

[7] C. L. Lawson and R. J. Hanson, Solving Least Squares Problems, Englewood 
Cliffs, New Jersey: Prentice-Hall, Inc., 1974. 


24 



[8] C. N. Guiar, “Using Ridge Regression in Systematic Error Corrections,” TDA 
Progress Report 42-92 , vol. October-December 1987, Jet Propulsion Laboratory, 
Pasadena, California, February 15, 1988. 

[9] G. H. Golub and C. F. Van Loan, Matrix Computations , Baltimore, Maryland: 
The Johns Hopkins University Press, 1989. 

[10] G. H. Golub, V. Klema, and G. W. Stewart, Rank Degeneracy and Least-Squares 
Problems , Technical Report TR-456, Department of Computer Science, Univer- 
sity of Maryland, College Park, Maryland, 1976. 

[11] W. H. Press, B. P. Flannery, S. A. Teukolsky, and W. T. Vetterling, Numerical 
Recipes , Cambridge, England: Cambridge University Press, 1986. 



Table 1. Systematic pointing error sources and model terms 


Error source 

Model function 

Cross-elevation error 

Elevation error 

Az Q collimation 

Pi c 

- 

Az encoder fixed offset 

P 2 cos (el) 

— 

Az/el skew 

P 3 sin (el) 

— 

Az axis tilt 

Pi sin (el) cos (az) 

— Pi sin (az) 

Az axis tilt 

P 5 sin (el) sin (az) 

P 5 cos (az) 

Source dec 6 

P 6 sin (az) 

Pe sin (el) cos (az) 

El c encoder fixed offset 

- 

P 7 

Gravitational flexure 

- 

P 8 cos (el) 

Residual refraction 

- 

P 9 cot (el) 

Az encoder scale error 

P 10 (az/360) cos (el) 

— 


Cross-declination error 

Declination error 

II A^ /dec axis skew 

— P\\ sin (dec) 

- 

HA axis tilt 

Pi 2 sin (HA) sin (dec) 

Pi 2 cos (1IA) 

HA axis tilt 

— Pi 3 cos (HA) sin (dec) 

P 13 sin (HA) 

HA feed offset 

-Pl4 

- 

Gravitational flexure 

Pi 5 cos (p)^ cos (el) 

-Pi 5 sin (p) cos (el) 

Declination feed offset 

- 

Pi 6 

Gravitational flexure 

Pit sin (p) cos (el) 

- 

Gravitational flexure 

- 

-Pi 8 cos (p) cos (el) 

Gravitational flexure 

-Pi 9 sin (el) 

- 

Gravitational flexure 

- 

P 20 sin (el) (el) 

HA encoder bias 

P 21 cos (dec) 

— 

rt Az refers to azimuth angle. 

d HA refers to hour angle. 


6 Dec refers to declination angle. 

c Uppercase P refers to parameter value. 

c El refers to elevation angle. 

f Lowercase p refers to paraleclic angle. 


Table 2. Applicable parameter sets to DSN 70-m and 34-m 
antennas 


Precision mode 


Computer 

command mode 

1 



l 

7 



9 

8 

0 



4 

11 



5 

12 



7 

13 



8 

14 



0 

21 

* 


- 


Table 3. Singular values for all-sky distributions 


Precision mode Computer command mod 


34.252 

34.252 

18.108 

18.096 

15.785 

11.814 

11.811 

11.807 

10.342 

10.312 

6.161 

5.641 

4.131 

2.623 

2.925 

0.819 

2.623 

- 

Condition number 

Condition number 

13.06 

4 1 .SO 
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Table 4. Computer command mode correlation matrix for all-sky distribution 


Parameter 

I 

2 

3 

1 

5 

7 

s 

9 

1 

1 .00 

- 0.07 

- 0.08 

0.00 

0.00 

0.00 

0.00 

0.00 

2 

— 0.07 

1 .00 

0.92 

0.00 

0.00 

0.00 

0.00 

0.00 

3 

- 0.08 

0.02 

1 .00 

0.00 

0.00 

0.00 

0.00 

0.00 

4 

0.00 

0.00 

0.00 

1 .00 

0.00 

0.00 

0.00 

0.00 

5 

0.00 

0.00 

0.00 

0.00 

1.00 

0.00 

0.00 

0.00 

7 

0.00 

0.00 

0.00 

a oo 

0.00 

] .00 

-0.81 

0.30 

S 

0.00 

0.00 

0.00 

0.00 

0.00 

- 0.81 

1 .00 

-0.75 

9 

0.00 

0.00 

0.00 

0.00 

0.00 

0.30 

- 0.75 

] .00 


Table 5. Singular values for the Voyager 2 data set Table 7. Theoretical standard deviations for the all-sky data 

— — __ set in millidegrees 

Precision mode Computer command mode — _ 





Parameter 

Precision inode 

Computer command mode 

57 . 22(5 

57.226 




J 1.530 

12.171 

1 

0.16 

0.82 

10.282 

10.282 

2 

- 

0.61 

5.912 

5.042 

3 

- 

0.69 

1.1 S 3 

1.483 

1 

- 

0.08 

0.334 

1.315 

5 

- 

o.os 

0.199 

0.199 

7 

0.17 

0.17 

0.0033 

0.0671 

8 

0.35 

0.35 

0.0028 

- 

9 

0.06 

0.06 

Condit ion number 

Condit ion number 

11 

0.23 


20,438 

853 

12 

0.08 

_ 


— 

13 

0.21 

_ 



14 

0.20 

_ 



21 

0.22 


Table 6. Theoretical standard deviations for the Voyager 2 data 

_ 



set in millidegrees 




Parameter Precision mode 

Computer command mode 

Table 8. Singular values for radio source distribution 


1 

2.19 

10.46 



2 


10.06 

Precision mode 

Computer command mode 

3 

- 

3.46 

15.951 

1 5 4 I '■'» 

4 

- 

0.20 

13.006 

13.157 

5 

— 

0.72 

9.885 

8.176 

8 

3.55 

3.55 

7.660 

7.504 

3.53 

3.53 

4.326 

1 238 

9 

11 

12 

0.03 

289.48 

0.72 

0.03 

3.550 

2.020 

2.603 

0.826 

13 



1.156 

0.283 

1 .89 

- 

O . S 25 


14 

252.03 

- 

Condition number 

Condition number 

21 

269.20 

" 

19.335 

54.600 
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Table 11. Singular values for reduced Voyager 2 parameter set 


Table 9. Theoretical standard deviations for the radio source 
data set in millidegrees 


Parameter 

Precision mode 

Computer command mode 

1 

0.33 

2.45 

2 

— 

1.68 

3 

_ 

1.96 

1 


0.13 

5 

— 

0.14 

7 

0.42 

0.42 

8 

1.11 

111 

9 

0.36 

0.35 

11 

0.44 

- 

12 

0.14 

— 

13 

0.44 

- 

11 

0.49 

- 

21 

0.57 



Table 10. Reduced parameter standard deviations for the 
Voyager 2 data set in millidegrees 


Parameter 

Precision mode 

Computer command mode 

1 

2.16 

- 

2 

_ 

0.29 

3 

- 

0.72 

4 

— 

0.20 

5 

_ 

0.72 

7 

3.55 

3.55 

8 

3.52 

3.53 

9 

0.03 

0.03 

12 

0.72 

- 

13 

1.S7 

- 

14 

1.30 

' 


Precision mode Computer command mode 


57.226 

57.226 

12.157 

10.282 

10.282 

8.893 

5.221 

4.845 

1.4S3 

1.483 

0.319 

1.315 

0.199 

0.199 

Condition number 

C on di t i on nu mber 

288 

288 


Table 12, Reduced parameter standard deviations for the 
Voyager 2 data set in millidegrees 


Parameter 

Precision mode 

Computer command mode 

2 

_ 

0.29 

3 


0.72 

4 

_ 

0.20 

5 

- 

0.53 

7 

- 

- 

8 

0.55 

0.55 

9 

0.03 

0.03 

12 

0.53 

- 

13 

0.19 

- 

14 

0.13 

— 


Table 13. Singular values for reduced Voyager 2 parameter set 


Precision mode 

Computer command mode 

56.955 

56.955 

S.952 

8.893 

S.137 

S.137 

4.900 

4.845 

1.335 

1 .335 

_ 

1.315 

Condition number 

Condition number 

43 

43 
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AZIMUTH, deg 



Fig. 2. Voyager 2 measurement set. 


Fig. 3. Radio source boresfght measurement set 









